function [c,c1] = show_basis_c1(dof_map_c0, dof_map_c1, V, T, weight, idx, d, d_plot)
% function [c,c1] = show_basis_c1(dof_map_c0, dof_map_c1, V, T, weight, idx, d, d_plot)

n_dof_c1 = max(max(dof_map_c1));
x = zeros(n_dof_c1,1);

if(idx > n_dof_c1)
    fprintf('too large index of d.o.f');
    idx = n_dof_c1;
end
x(idx) = 1;

n_dof_c0 = max(max(dof_map_c0));
c = zeros(n_dof_c0,1);
c1 = zeros(21,size(T,1));
for k = 1:size(T,1)
    c_loc = weight(:,:,k)*x(dof_map_c1(:,k)); 
    c(dof_map_c0(:,k)) = c_loc; 
    c1(:,k) = c_loc;
end

[vtx, tri, u] = heval(dof_map_c0, V, T, c, d, d_plot);

trisurf(tri, vtx(:,1), vtx(:,2), u);